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Abstract 

Quantum superintegrable systems in two dimensions are obtained from their 
classical counterparts, the quantum integrals of motion being obtained from 
the corresponding classical integrals by a symmetrization procedure. For each 
quantum superintegrable system a deformed oscillator algebra, characterized by 
a structure function specific for each system, is constructed, the generators of 
the algebra being functions of the quantum integrals of motion. The energy 
eigenvalues corresponding to a state with finite dimensional degeneracy can 
then be obtained in an economical way from solving a system of two equations 
satisfied by the structure function, the results being in agreement to the ones 
obtained from the solution of the relevant Schrodinger equation. Applications 
to the harmonic oscillator in a flat space and in a curved space with constant 
curvature, the Kepler problem in a flat or curved space, the Fokas-Lagerstrom 
potential, the Smorodinsky-Winternitz potential, and the Holt potential are 
given. The method shows how quantum algebraic techniques can simplify the 
study of quantum superintegrable systems, especially in higher dimensions. 



1 Introduction 



The idea of studying the properties of physical systems exhibiting degenerate 
energy levels through use of their symmetries has been exploited since the early 
days of quantum mechanics. In these cases the symmetry algebra Q of the sys- 
tem has to be determined, which is a finite-dimensional Lie algebra containing 
ladder operators which connect all the eigenstates with a given energy, while the 
Hamiltonian of the system is related to the Casimir operators of the algebra. 
The set of eigenstates with given energy provides a basis for an irreducible rep- 
resentation (irrep) of the algebra. The energy eigenvalues are then determined 
by the eigenvalues of the Casimir operators of the algebra in the corresponding 
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irreps. The determination of the symmetry algebra of a given hamiltonian is a 
quite difficult task, which is not always dictated by an a priori obvious proce- 
dure of searching for it. Examples of physical systems having known symmetry 
algebras are the N-dimensional isotropic harmonic oscillator and the Kepler 
(Coulomb) system, bearing the symmetries SU(N) and SO(N+l) respectively. 

The classical counterparts of the quantum isotropic harmonic oscillator and 
the quantum Kepler problem have another interesting property. They are max- 
imally superintegrable systems in N dimensions. Superintegrable systems in N 
dimensions have more than N independent classical constants (also called in- 
tegrals or invariants) of motion, while maximally superintegrable systems in N 
dimensions have 2N— 1 independent classical constants of motion, N of which 
are integrals in involution (the Poisson bracket of each pair of them is zero). 
This property of the maximally superintegrable systems implies, in classical 
mechanics, that every closed and bounded trajectory is a periodic trajectory. A 
detailed review of superintegrable systems in 2 dimensions is given in H , while 
examples of superintegrable systems in 3 dimensions can be found in |J . 

Higgs H and Leemon || have shown that in the case of a N dimensional 
system moving in a space with constant curvature the isotropic harmonic os- 
cillator and the Kepler problem are still maximally superintegrable systems, 
both in classical and quantum mechanics. Furthermore, they have shown that 
the quantum counterparts of these systems can be described by symmetry al- 
gebras isomorphic to SU(N) and SO(N+l) respectively. Additional examples 
of superintegrable classical systems are the Fokas-Lagerstrom potential |7| , the 
Smorodinsky-Winternitz potentia l jR ||, [n], the Holt potential Jl2f , the 
Hartmann potentials [ p"3"[ |l4| |l5| , |16| . The problem of quantum integrability 
and its connection to classical integrability is currently under active investiga- 
tion 0,0,0,0. 

In several of the above mentioned cases (see |5|, |6[ , for example) the property 
of the classical and quantum superintegrability of a physical system coincides 
with the existence of a symmetry algebra of the system. The energy levels of 
the system can then be determined by purely algebraic means. However, the 
identification of a symmetry algebra is not always easy. Furthermore, in some 
cases (see jl6) , for example) it seems that the usual Lie algebras do not suffice for 
this purpose. The recent introduction of quantum algebras |2l], ^2, 23, 24[ (also 
called quantum groups) opens new possibilities in this direction. Quantum alge- 
bras are nonlinear deformations of the usual Lie algebras, to which they reduce 
when the deformation parameter q goes to 1. They were initiated as a mathe- 
matical tool issued from the study of the quantum inverse problem, the Yang 
Baxter equation and conformal field theories (see for a collection of origi- 
nal papers). There are already some indications that quantum algebras might 
be useful as symmetry algebras of certain superintegrable systems. The Higgs 
algebra (i.e. the symmetry algebra of the Kepler problem in a 2-dimensional 
space with constant curvature studied in j^]) can be approximated to second 
order by the quantum algebra SU g (2) j^6|. The symmetry algebra of the Hart- 
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mann potential, for which usual Lie algebras seemed insufficient |fq ], has been 
identified as the quadratic Hahn algebra QH(3) ^7|. The quadratic Hahn al- 
gebra QH(3) has also been found to describe the symmetry of the anisotropic 
singular oscillator (a 3-dimensional harmonic oscillator with an additional term 
~ f/(r 2 sin 2 6*) The quadratic Hahn algebra QH(3) is a special case of 

the general quadratic Askey- Wilson algebra QAW(3), which is the dynamical 
symmetry of the potentials having eigcnfunctions described by classical polyno- 
mials Letourneau and Vinet Q have constructed the quadratic algebra 
describing the case of a harmonic oscillator potential with a 2 to f frequency 
ratio. Another example of a system described by a nonlinear algebra is the 
generalized Kepler (Coulomb) system From these examples it becomes 
clear that nonlinear algebras can be useful in the description of integrable and 
superintegrable systems. 

In this paper we focus attention on the simplest superintegrable systems, the 
two dimensional ones, and we propose a method of determining their dynamical 
symmetries and calculating their spectra by purely algebraic means. It turns out 
that several quantum superintegrable systems in 2 dimensions can be described 
in terms of appropriate generalized deformed oscillators, which allow for the 
direct determination of the energy levels and their degeneracies without any 
need of solving the Schrodinger equation. A preliminary study of the proposed 
method can be found in ref. |52|, where the method has been used in two 
cases of potentials, the symmetric harmonic oscillator in a curved space and the 
asymmetric oscillator with a 2 to 1 frequency ratio. In this paper we study most 
of the known 2 dimensional superintegrable cases. 

It is known that the classical (non-deformed) algebras can be constructed 
using the harmonic oscillator algebra {a, a + , N} as the basic underlying struc- 
ture. For quantum algebras, and nonlinear algebras in general, deformed os- 
cillators have to be used. Biedenharn [[33) and Macfarlane [[m) constructed the 
q-deformed oscillator appropriate for the Schwinger realization of the quantum 
algebra SU g (2). Many other deformed oscillators can be found in the literature. 
We mention the Q-oscillator introduced by Arik and Coon (3f| and Kuryshkin 
[ B6) , the two-parameter deformed oscillator |}7| [|q], the parafermionic oscilla- 
tor (39) and its q-deformation fyj, the parabosonic oscillator []39| a nd its 
q-deformation [[lo[ [ft), the generalized q-deformed fermionic algebra [|42|. (The 
q-deformed fermionic algebra jl3| has been proved to be equivalent to the usual 
fermionic algebra ]44[ .) The common feature of all these deformations is their 
structural similarity. In all cases an appropriate Fock basis can be constructed, 
leading to a matrix representation of the algebra. 

The structural similarity of the various deformed oscillators implies that all of 
them can be described in a unified framework. Among the various alternatives 
we mention here the generalized deformed oscillator ^5|, the Odaka-Kishi- 
Kamefuchi unification scheme ]4l|] , the pioneering work of Jannussis et al. on 
the bozonization method and the generalized Q-deformed oscillator pTy , the 
Beckers-Debergh unification scheme [Eq, and the Fibonacci oscillator scheme 
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p9| , while a general treatment of qualgebras is given recently by Fairlie and 
Nuyts@. 

Among the various equivalent descriptions of the deformed oscillators, in this 
paper we use the deformed oscillator algebra p5[, already used for the study of 
the energy spectra of one dimensional systems [pT|, [d2| [d3[ . 

In section 2 of the present paper we consider the classical superintegrable 
systems in 2 dimensions, for which we show that the relevant Poisson bracket 
induced algebra has a structure similar to the deformed oscillator algebra. In 
section 3 a working hypothesis for the quantum superintegrable systems in 2 
dimensions is proposed, which leads to the calculation of the energy eigenvalues 
and their degeneracies by purely algebraic means. This hypothesis is applied to 
several quantum superintegrable examples in section 4, while section 5 contains 
discussion of the present results and plans for further work. 

2 Classical superintegrable systems in two di- 
mensions 

Consider a classical system with two degrees of freedom, described by the Hamil- 
tonian: 

H = H(x,y,p x ,p y ). (1) 

If the system is superintegrable there are two independent additional integrals 
of motion / and C, such that: 

{H,I} PB ={H,C} PB = 0, and {/, C} PB = F(H, I, C), (2) 

where { , } PB denotes the Poisson bracket and F = F(H, I, C) is a constant of 
motion which depends on the three independent constants of motion H, I, C. A 
superintegrable system in two dimensions is necessarily a maximally superinte- 
grable system, which means that all finite classical trajectories are closed and 
periodic. Integrable and superintegrable systems in 2 dimensions have been re- 
viewed in ^ , while in j^] a systematic study is given of superintegrable systems 
in 3 dimensions which possess invariants that are quadratic polynomials of the 
canonical momenta. 

Maximally superintegrable systems possess, by definition, the maximum 
number of independent classical invariants. Therefore any other integral can 
be expressed as a function of the basic integrals H, I, C. As a result we can in 
general choose two new integrals of motion: 

L = L(H,I,C), and A = A(H, I, C), 

such that: 

{L,A} PB =B, {L,B} PB = -A. (3) 
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After a calculation we can prove that: 

B 2 + A 2 = G(H, L), 

where G{H, L) is some function depending only on the integrals of motion iJ, 
L, and 

i f)(~< 

{ A,B} PB = <P(H,L) = -- — . (4) 

The structure of the algebra defined by eqs ([||4|) has many similarities to 
the algebraic structure of the deformed oscillator given in references |ff5 
where L is some kind of number operator, while A, B are like the creation and 
annihilation operators. The deformed oscillator algebra is a non abelian algebra. 
Therefore it is quite natural to attempt studying the quantum superintegrable 
systems by applying some similar procedure in order to calculate their quantum 
properties (eigenvalues and eigenvectors), since the corresponding properties of 
the deformed oscillators have been already studied. 



3 Quantum superintegrable systems in two di- 
mensions 

The question of intcgrability in quantum mechanics is under investigation. Many 
authors |l7], [l8|, [H], have investigated several aspects of the extension of 
integrability from classical to quantum mechanics. We should recall at this point 
that each quantum system with discrete energy spectrum can be considered as 
a quantum integrable system fTi| . 

In this paper we consider a two dimensional quantum system described by a 
hamiltonian H acting on a Hilbert space TL. The hamiltonian is an autoadjoint 
operator with range dense in the space TL. All the operators defined in this 
section are supposed to be generated by nonlinear combinations of the basic 
algebra of generators x,p Xl y,p v satisfying the usual commutation relations: 

= [y,p y ] = i, other commutators = 0. 

If the system is integrable, then there is a second autoadjoint operator I 
commuting with the hamiltonian and having a range dense in TL 

[H,I]=0, (5) 

the operators H and / being linearly independent. The commutativity of these 
operators implies that there is a family of common eigenvectors for both op- 
erators. Let us label the common eigenvectors of these operators using their 
corresponding eigenvalues, thus obtaining a family of vectors in the Hilbert 
space TL. Let us further suppose that this family spans the whole Hilbert space 
and that there is no simultaneous degeneracy for both labels of the common 
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eigenfunctions (i.e. there could be degeneracy in each label separately, but not 
in both labels simultaneously). This assumption defines some kind of complete- 
ness of the system of the two quantum integrals H, I in involution. The present 
assumption is also consistent with the fact that the ranges of the operators H 
and I are supposed to be dense in the ambient Hilbert space 7i. 

Weigert |ljj has proved that there is another quantum integral I' possessing 
the same properties as /. This is permissible in quantum mechanics but not 
in classical mechanics. In this sense the choice of complete operators H, I' 
instead of H, I means a different way of labelling of the base spanning the 
Hilbert space. Here we shall treat the part of the Hilbert space corresponding 
to a discrete spectrum. In some cases the whole of the Hilbert space can be 
described by a discrete basis (as in the case of the harmonic oscillator), while 
in other cases (as in the case of the Coulomb (Kepler) potential) a part of the 
Hilbert space is described by a discrete basis (the space of energy eigenvectors 
with negative energy eigenvalues). 

The system is called superintegrable, by analogy to the classical definitions, 
if there is a third operator C, linearly independent from H and /, with range 
dense in the Hilbert space and commuting with H but not commuting with / 

[H,C]=0, [J,C]^0. 

In this paper we propose the following working hypothesis: 
Hypothesis: Let us consider the superintegrable systems for which we can 
construct an associative algebra: 



(6) 



N = 


M{H,I,C), 


N+ = 


M, 


A = 


A(H,I,C), 


W,A = 


-A, 


A+A = 


$(H,A0, 




0, 


is a real positive function 


definite for x 


<f>(E,0) 


= 0. 



(<) 



From the above equations we can prove that: 

[N,A+]= A+, 

AA + = $(H,Af+l), 

If this construction is possible we can then define the Fock space for each energy 
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eigenvalue: 

H\E,n>= E\E,n>, 

Af\E,n> = n\E,n >, n = 0, 1 . . . , 

A\E,0>= 0, 

\E,n>= (^^=)j(A + ) n \E,0>, 

where 

[0]! = 1, [n]l = $(E,n)[n-l]L 

In the case of the discrete energy eigenvalues, for every energy eigenvalue E 
there is some degeneracy of dimension + 1 . Therefore the dimensionality of 
the Fock space corresponding to that energy eigenfunction should be equal to 
N c i + 1. This is equivalent to the condition: 

$(E,N d + l) = 0. (8) 

As we shall see in the examples given in the following section, the two conditions 
([J) and fig), and the positiveness of the structure function $(E, x) suffice in order 
to determine the energy spectrum of the quantum maximally superintegrable 
systems. 

There are only a few quantum 2-dimensional superintegrable systems known. 
All the examples studied in this paper have a classical counterpart. 

4 Examples of quantum superintegrable 
two-dimensional systems 

In this section we shall apply the hypothesis of the previous section in or- 
der to determine the energy spectrum of some quantum superintegrable two- 
dimensional systems using purely algebraic methods. 

4.1 Harmonic oscillator in a flat space 

The two-dimensional symmetric harmonic oscillator in euclidean coordinates is 
described by the hamiltonian: 

H= 1 -(pl+pl) + l^(x 2 +y>). (9) 

The following Fradkin operators H| can be defined: 

B = S XX - S yy = (pi + uj 2 x 2 ) - (p 2 v + u 2 y 2 ) , S xy = p x p y + cu 2 xy. (10) 

The operator B in (|Io| ) is the quantum mechanical analogue of a third con- 
stant of motion in the sense of classical hamiltonian mechanics , the second 
one being the angular momentum operator: 

L = xp y -yp x , (11) 
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since: 

[H,L] = [H,B]=0, (12) 

while the operators B, L do not commute, but they form a closed algebra with 
the operator S xy : 

[L,B}=4iS xy , [L,S xy ] = -iB. (13) 

The above relations suggest the possibility of expressing the two-dimensional 
harmonic oscillator algebra by using the deformed oscillator formulation: 

A+ = ^ + iS xy , (14) 

A=--iS 

where u is a constant to be determined and 

W,A+]= A+, 

{Af,A} = -A, 

A+A= H 2 -u 2 {L-l) 2 (15) 

= H 2 -uj 2 {2N+2u- l) 2 

= *{H,Jf), 

where the function <fr(E,x) is given by: 

$(E, x) = E 2 - lu 2 (2x + 2u - l) 2 , (16) 

and we can see that 

AA + = &(H,SS+1). 

The existence of a finite dimensional representation of the oscillator algebra 
is equivalent to the existence of a maximum number N + 1 which is a root of 
the structure function, with N being the dimensionality of the algebra represen- 
tation, coinciding with the dimensionality of the appropriate Fock space. This 
restriction, combined with the annihilation of the structure function for x = 0, 
is written as: 

*(E,0) = 0, 
$(E,N + 1) = 0, (17) 
$(E,x)>0 for x = 1,2,..., N. 

Solving this system of two equations with two unknowns, E and u, one obtains 
the eigenvalues of the harmonic oscillator in a flat space: 

N 

u=- — , E = E N = lo{N+ 1). (18) 
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The angular momentum values allowed for each energy level can then be deter- 
mined by inserting the value of the constant u just obtained into the first of eq. 
(p~4|) , the result being: 

L = -N,-N + 2,...,N-2,N. (19) 

The structure function of the deformed oscillator is calculated to be: 

$(E Ni x) = 4oj 2 x(N + 1 - x). 

Clearly the above energy ( [L8| ) and angular momentum spectra (|l9|) are the 
same as the ones obtained by classical means, i.e. by solving the appropriate 
Schrodinger equation, or by using the SU(2) symmetry. The existence of a 
finite dimensional algebra representation should be attributed to the existence 
of stable periodical trajectories in the corresponding classical case. 



4.2 Harmonic oscillator in a space with constant curvature 

Higgs H has studied the symmetries of a harmonic oscillator in a non-flat space, 
a space with constant curvature in particular. A typical example of such a space 
is the surface of the sphere in a three dimensional space. 

The curved space is geometrically described by the metric: 



dx 2 + dy 2 + X(xdy — ydxY 
{l + X{x 2 +y 2 )f 



ds 2 = 



the flat space corresponding to A = 0. The harmonic oscillator in this space is 
defined in ref. J5| by the Hamiltonian: 

H= 1 -(^ + n 2 + XL 2 ) + ^-(x 2 + y 2 ), (20) 



where the angular momentum operator L is given by eq. ( |11| ) and 

TTx = Px + | {x (xp x + yp y ) + (xp x + yp v ) x) , 
% = Py + 2 (V ( x Px + VPy) + ( x Px + VPy) v)- 

By analogy to the harmonic oscillator in a flat space, Higgs [|| has defined the 
Fradkin-like operators: 

B = S XX - Syy = {tt 2 x + uj 2 x 2 ) - (n 2 y + uj 2 y 2 ) , (22) 

S xy = -^{^x,^y} + u 2 xy, (23) 

which are symmetrized versions of eq. (|To|). 

The commutators of H with L and B are given in eq. (|l2|), while the 
operators B, L do not commute and their commutator is the same as in equation 



(^3[), where the coordinates of the momentum p should be replaced by the 
coordinates of the extended momentum w, defined by eq. and already used 
in the symmetrized Fradkin operators given above. 

The operators H, L, B, S xy define again a closed non-linear algebra as in the 
flat harmonic oscillator case. Therefore we have another example of a superin- 
tegrable system in a non-flat space. 

In the present case we can also define the corresponding deformed oscillator 
as in eq. (^J). The deformed algebra is then completely defined by the structure 
function: 

*(E,x) = 

E 2 -(uj 2 + ^ + XE^j (2x + 2u - l) 2 + £ (2x + 2u- l) 4 . (24) 

Following the same methodology as in the case of the harmonic oscillator 
in a flat space, we can generate the corresponding Fock space for the curved 
harmonic oscillator. By assuming then the existence of a finite dimensional 
deformed algebra representation the restrictions corresponding to eq. ([[t]) are 
valid. These equations determine the energy eigenvalues: 



E = E n = ] Jlu 2 + ^(N+1) + ^(N + 1) 2 , (25) 

while the constant u turns out to have the same values as in eq. (|l^). The 
angular momentum eigenvalues are again given by eq. (|l9|). 

Another interesting point arises from the comparison between the structure 
function ( |lrj| ) in a flat space and eq. (|24|) in a curved space: The geometry of 
the space affects the algebra characterizing the harmonic oscillator. For A = 
eq. (|i|) reduces to eq. (|l6|), as it should. 

Finally, the structure function can be written as 

$(E N ,x) = 4x(N + l-x) (\(N+l-x) + yju 2 + A 2 /4) (\x + y/u 2 + A 2 /**) . 

The symmetries of the harmonic oscillator in a curved space have been studied 
in ref. |p5| using the notion of the quadratic Racah algebras QR(3). 



4.3 The Kepler problem in a curved space 

The case of the Kepler problem in a space with constant curvature has been 
studied by Higgs || . The hamiltonian is given by: 

H = l(n 2 x + n 2 y + XL 2 )-^ r = v^+^, (26) 

where the angular momentum operator L is given by eq. Jll] ) and the ir x , tt v 
are defined in eq. (|2l|). 



10 



The Runge-Lenz vectors in the curved space can be defined by: 



R, =-~{L,TTy} + tl~, Ry=^{L,TT X } + fl^. (27) 



This system is a quantum superintergrable system in a curved space because: 

[H,L]=0, [H,R x } = 0, 
and the operators L, R x , R y form a closed algebra: 

[L, R x ] = iR y , [L, R y ] = —iR x . 
Using the same hypothesis as previously we can define the deformed oscillator 



algebra: 




N = 


L — u, 


A+ = 


Rx + iRy , 


A = 


Rx ^Ry i 


A+A = 


H 2 + 2H (L- 1/2) 2 




2H{N + u- 1/2) 2 







The structure function in this case is defined by: 

$(E,x) = + 2E(x + u-l/2) 2 

—A (x + u — 1/2) 2 ((.x +u- 1/2) 2 - 1/4) 
The solution of eqs (p7|) is given by: 



(28) 



N 2u 2 N(N + 2) , . 

The permitted eigenvalues of the angular momentum operator L are given by: 

N N N N 

This means that the symmetries of the Kepler problem are compatible with the 
existence of angular momenta equal to 0, 1/2, 1, 3/2,. ... In physical situations, 
however, only integer angular momenta appear, which means that N — 2n. In 
this case the spectrum given by eq. (|9|) is the same to that obtained by Higgs 
[|J. In the case of zero curvature, i.e. A = 0, we obtain the usual Coulomb 
energy spectrum. 
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The structure function corresponding to the Kepler problem is given by: 



, 4u 2 JN + 1- 2x)' 
®{E N ,x)=x(N + l-x) [ tAT " +A- V 



(N+l) 2 4 
Zhedanov [E6| has proven that the nonlinear algebra of eq. (E8h can be approx- 



imated to second order by the SU g (2) algebra |33 



The symmetries of the Kepler potential in a curved space have been studied 
in ref. []56| using the notion of the quadratic Racah algebras QR(3). 

4.4 Fokas-Lagerstrom potential 

In classical mechanics the superintegrable system described by the Hamiltonian: 



has been studied by Fokas and Lagerstrom |]]. This system has two additional 
classical invariants of motion, 

3 2 

J = pl+x 2 1 and C = (xp y - yp x )p 2 y + - , (31) 

the second of which (C) is a cubic function of the coordinates. The quantum ver- 
sion of the hamiltonian ([so]) corresponds to a quantum superintegrable system 
with two additional integrals: 

1 y 3 p \xy 2 PyJ 
J = pl + x 2 , and B = - {xp v - yp x ,p 2 y } + -— ^ fT"^' 

where { , } is the usual anticommutator. It is clear that the quantum integral B 
is the symmetrized version of the classical integral C . From the above definitions 
we can verify that: 

[H,J] = 0, [H,B]=0, 
[J,B] = R, [J,R]=4B, 

and 

56 qo 

[R, B] = 8 J 3 - 36 J 2 H + 48JH 2 - 16H 3 + — J — —H, 

9 9 

R 2 - 4B 2 = 4J 4 - 2AJ 3 H + A8J 2 H 2 - 32JH 3 + 



+ mj2_§m JH + A SH 2 + f. 
From the above closed algebra we can define: 

M = J/2 -u, A + = B + R/2, A = B~ R/2, 
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where u is a constant to be determined. These operators correspond to a de- 
formed oscillator algebra: 



[Af,A+} = A+, [Af,A] = -A, 

A+A = -(J -1){2H - J+ 1)(6H - 3J + 1)(6H - 3J + 5) 
9 

= i(2AT-l + 2u)(2£r-2«+l-2A0 (32) 
(6H - 6u + 1 - 67V)(6F - 6u + 5 - 67V) 

A4+= $(i2",7V+l). 

The corresponding structure function is defined by: 

$(E,x)= ±(2x-l + 2u)(2E-2u+l-2x) 

{GE-6u+l- 6x)(6E - 6u + 5 - 6x). 

The existence of a finite representation of the algebra for each energy eigen- 
value implies that the structure function satisfies eq. ([l7]) and it is a positive 
function. Therefore we can find the possible energy eigenvalues having degen- 
eracy equal to N + 1: 

Case i) 

it = 1/2, and E N =N + 1, 
corresponding to the structure function: 

${E N , x) = 16x(7V + 1 - x) \N + | - xj \ N + | - x 

Case ii) 

tt=l/2, and E N =N + 2/3, 
corresponding to the structure function: 



&(E N , x) = \Qx{N + l-x){N+--xj [N + 
Case iii) 

u = 1/2, and E N = N + 4/3, 
corresponding to the structure function: 



®{E N , x) = \Qx{N + l-x)[N+--x) [N + 
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In all cases the degeneracy is determined by J = 2(JV + u). Since M obtains 
the N + 1 values 0, 1, . . . , N, as a result J also obtains N + 1 values. 

The Hamiltonian (30) corresponds to the linear combination of two harmonic 
oscillators which have the above described energy spectrum. This case has a 
special significance, since it is an example of a superintegrable potential which is 
not a separable one in two different coordinate systems. The proposed method 
does not depend on the separability of the variables in two systems. 



4.5 Smorodinsky— Winternitz potential 

The classical superintegrable Smorodinsky-Wintcrnitz system 
^. pi |ll) corresponds to the Hamiltonian: 

H=\{p 2 x +p 2 y )+k{x 2 +y 2 ) + ^. (33) 
This system has two additional classical invariants of motion, 

v 2 

T = p 2 y + 2ky 2 , and C = x 2 p 2 y + y 2 p 2 x - 2xyp xPv + 2c^, (34) 

the second of which (C) is a quartic function of the coordinates. Evans |ll[] has 
proved that the Winternitz-Smorodinsky potential in N dimensions is an exam- 
ple of a superintegrable system. The quantum version of the hamiltonian ( p3[ ) 
corresponds to a quantum superintegrable system with two additional integrals: 

T = p 2 y + Iky 2 , and 5 = x 2 p 2 y + y 2 p 2 x - {xy , p x p y } + 2c^ ■ 

It is clear that the quantum integral B is the symmetrized version of the classical 
integral C. From the above definitions we can verify that: 

[H,T} = 0, [H,B] = 0, 

and 

[T, B] = R, [T, R] = 32kB + 8T 2 - \QHT - 16fc, 

[R, B] = IQBT ~ 16BH + 32(c - 1)T + 8R + 32H, 

R 2 = 32kB 2 + 224kB + 32(c - 1)T 2 + 64HT + IQRT - IQRH - 48# 2 
+ 1&BT 2 - 32BTH + 192fc(c - 1). 

From the above closed non-linear algebra we can define: 



TV = 


J—T 


- u, 








A+ = 


AkB+ y 




- T 2 


-2HT- 


2k, 


A = 


AkB — y 




- T 2 


-2HT- 


2k, 
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(35) 



where u is a constant to be determined. These operators correspond to a de- 
formed oscillator algebra: 

W,A+]= A+, [tf,A]=-A, 

A + A= 24H 2 k + 3 • 22 Hk§ + 36fc 2 

-96 c fc 2 - 2i iJ 2 v^T - 88 H kT 
-3-2ifc§T + 2^ cklT + AH 2 T 2 
+3 • 2i Hy/kT 2 + UkT 2 ~ 16 ckT 2 
-4HT 3 - 23 VkT 3 + T 4 
= $(H,Af), 
AA+ = ${H,N + 1). 

The corresponding structure function can be factorized as: 

$(E, x) = 1024fc 2 (x - (u + §)) far - (t*+ 3)) 

The existence of a finite representation of the algebra for each energy eigenvalue 
implies that the structure function satisfies eq. (|l7]). The positiveness of the 
structure function for every < x < N implies: 

3 

U = —^ 

while the energy eigenvalues are given by: 



with — i < c. The structure function is given by: 



3>(E N ,x) = 

= I024k 2 x (x + |) (N + 1 - x) (N + 1 + 
If the following restriction is valid: 



yT+gc 
2 



8 

the following energy eigenvalues are also permitted: 



-l<c<l, (36) 



corresponding to the structure function: 



®{E N , x) = 1024fc 2 x ( x + ^]{N + 1 - x) (n + 1 - ^ j~ 8c 
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In both cases the degeneracy of the levels is determined by T — \/32k(Af—u). 
Since Af obtains the N + 1 values 0, 1, . . . , N, as a result T also obtains N + 1 
values. 

It is worth noticing that in the case of solving the problem using the Schro- 
dinger equation, the restriction (^) is introduced by the assumption that the 
eigenfunctions should be square integrable functions on the plane (x,y). In 
the Schrodinger equation solution the additional restriction of finiteness of the 
potential energy restricts the choice of c to positive values only. 



4.6 The Holt potential 

The classical superintegrable Holt |L2| system corresponds to the Hamiltonian: 

H = \{&+1® + + (37) 

This potential is a generalization of the harmonic oscillator potential with a 
ratio of frequencies 2:1. This system has two additional classical invariants of 
motion, 



2(5 

T = p 2 +8y 2 1 and C = p 2 x p y + 8xyp x - 2x 2 p y + — p y , (38) 



the second of them (C) being a cubic function of the momenta. The quan- 
tum version of the hamiltonian (37) corresponds to a quantum superintegrable 
system with two additional integrals: 

26 

T = p 2 + 8y 2 , and B = p 2 p y + 4 {xy,p x } - 2x 2 p y + —p y . 

It is clear that the quantum integral B is the symmetrized version of the classical 
integral C . From the above definitions we can verify that: 

[H,T} = 0, [H,B]=0, 

and 

[T, B] = R, [T, R] = 32B, 
[R, B] = -96 + 2565 - 64H 2 + 128HT - 48T 2 , 
R 2 - 32B 2 = 102477 - 704T + 512<5T - 128TH 2 + 128T 2 H - 32T 3 . 
From the above closed non-linear algebra we can define: 

TV = -= - u, A + = 8B + V2R, A = 8B- V2R, 
v/32 
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where u is a constant to be determined. These operators correspond to a de- 
formed oscillator algebra: 

W,A+]= A+, [N,A] = -A, 
A+A = 2 6 (T — 2a/2) 



H-Z + y/2 + Jm\lH-Z + ,/2-Jm) (39) 



AA+ = $(H,tf+l). 
The corresponding structure function is defined by: 

9(E,x)= 2f((x + u)-|) 

_(* + „) + § + 

The existence of a finite representation of the algebra for each energy eigenvalue 
implies that the structure function satisfies eq. (|l7j). Therefore we can find the 
possible energy eigenvalues having degeneracy equal to N + 1 : 



I, and E N = V8 (n + 1 X ' 



T \ 4 

where (1 + 85) > 0. The corresponding structure function is: 

VI + 8<5 



x) = 2^x(N + l-a;)lJV + l- a;+ ^ 
In the special case where — | < <5 < I there are energy eigenvalues given by: 



u = ~ and Bjv = V8 TiV + 1 - ^— ±J 
and the structure function is: 



$(£jv, x) = 2-x(7V + 1 - x) \ N + 1 - 

which is positive for < x < N if — | < <5 < |. 

In both cases the degeneracy of the levels is determined by T = V32(A/" + u). 
Since N is obtaining the N + 1 values 0, 1, . . . , iV, as a result T also obtains 
N + 1 values. The quantum Holt potential has also been studied recently by 
using quadratic algebras by Letourneau and Vinet pQj. 
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5 Discussion 



In this paper, starting from classical superintegrable systems, we have shown 
that the corresponding quantum systems are superintegrable ones, the quantum 
integrals (quantum constants of motion) being obtained from the classical ones 
using a symmetrization procedure. Furthermore, the quantum superintegrable 
systems can be described in terms of a deformed oscillator algebra. The op- 
erators of the deformed oscillator algebra are constructed from the quantum 
integrals. The deformed oscillator algebra is characterized by a structure func- 
tion &(E,N), which takes a specific form for each superintegrable system. The 
eigenvalues of the energy and their degeneracies are determined in an economical 
way directly from equations satisfied by the structure function, the results being 
in agreement with these coming from the independent solution of the relevant 
Schrodinger equation. 

A few comments and some open problems are now in place: 

i) In all of the examples considered in this paper, quantum superintegrability 
is induced by classical superintegrability, the quantum integrals of motion being 
symmetrized versions of the corresponding classical integrals. The extend to 
which classical superintegrability implies in general quantum superintegrability 
as well, has to be tested. 

ii) In all of the examples considered in this paper, quantum superintegrability 
manifests itself in the degeneracy of the energy levels, a fact already noticed Q. 

iii) The hypothesis that to each superintegrable system corresponds a de- 
formed oscillator algebra, i.e. 

superintegrability — > deformed oscillator algebra 

is an exact proposition in the classical case, as shown in section 2. In the 
quantum case, however, a general formal proof is still lacking. In section 3 
a working hypothesis was made, which was proved successful in the examples 
considered in section 4. 

iv) The list of two-dimensional quantum superintegrable systems given in 
this paper is not exhaustive. There are classical superintegrable systems for 
which the quantum superintegrability has to be proven, as, for example, the 
Calogero system ([Q, eq. (3.5.9)), which possesses a sixth order invariant. Fur- 
thermore, there are two-dimensional systems for which the quantum superinte- 
grability has been shown, but the determination of the corresponding deformed 
oscillator algebra requires heavy computation, as, for example, the Winternitz- 
Smorodinsky potential of ref. |9J, given also in @, eq. (3.2.36). 

v) The extension of the present method to three-dimensional quantum super- 
integrable systems is under investigation. It should be mentioned that classical 
superintegrable systems in 3 dimensions having invariants which are quadratic 
polynomials in the canonical momenta have been recently studied in 

vi) Another interesting point is the semiclassical study of two-dimensional 
superintegrable systems. The algebra characterizing the two-dimensional classi- 
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cal superintegrable systems, studied in section 2, can be quantized by using the 
correspondence { Poisson bracket — > commutator }. Through this procedure, 
from the classical algebra a quantum deformed algebra is obtained, which is the 
semiclassical counterpart of the exact quantum oscillator algebra considered in 
this paper, these two algebras being slightly different. 

vii) The example of the Fokas Lagerstrom potential, which is the oscilla- 
tor with ratio of frequencies 1:3, shows that quantum superintegrability implies 
a dynamical symmetry. This example was examined using algebraic methods 
for the first time. All the other examples have been already studied by other 
authors, as it has been indicated in the text. The difference of the proposed 
treatment is that we do not use the separability in two coordinate systems in 
order to calculate the dynamical symmetries. The set of two-dimensional sys- 
tems separable in two different coordinate systems is a subset of the class of the 
superintegrable systems in two dimensions. The study of other superintegrable 
systems non-separable in more than one coordinate system seems to be very 
interesting. A class of such systems already known consists of the oscillators 
with rational ratio of frequencies. 

viii) Many of these examples f3Q, B9, p6j have been studied by using cases of 



quadratic Askey- Wilson algebras QAW(3) |29|. An open problem is if the gen- 
eral quadratic Askey- Wilson algebra can be expressed by a deformed oscillator. 
One can also notice that the algebra of the Fokas-Lagerstrom problem is a cu- 
bic algebra, while all the other examples in this paper correspond to quadratic 
algebras. 
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